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Abstract 

We consider Horava gravity within the framework of the effective field theory (EFT) of dark energy and modified 
gravity. We work out a complete mapping of the theory into the EFT language for an action including all the 
operators which are relevant for linear perturbations with up to sixth order spatial derivatives. We then employ 
an updated version of the EFTCAMB/EFTCosmoMC package to study the cosmology of the low-energy limit of 
Horava gravity and place constraints on its parameters using several cosmological data sets. In particular we use 
cosmic microwave background (CMB) temperature-temperature and lensing power spectra by Planck 2013, WMAP 
low-f! polarization spectra, WiggleZ galaxy power spectrum, local Hubble measurements, Supernovae data from 
SNLS, SDSS and HST and the baryon acoustic oscillations measurements from BOSS, SDSS and 6dFGS. We get 
improved upper bounds, with respect to those from Big Bang Nucleosynthesis, on the deviation of the cosmological 
gravitational constant from the local Newtonian one. At the level of the background phenomenology, we find a 
relevant rescaling of the Hubble rate at all epoch, which has a strong impact on the cosmological observables; at the 
level of perturbations, we discuss in details all the relevant effects on the observables and find that in general the 
quasi-static approximation is not safe to describe the evolution of perturbations. Overall we find that the effects of 
the modifications induced by the low-energy Horava gravity action are quite dramatic and current data place tight 
bounds on the theory parameters. 
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1 Introduction 

In their quest to find a quantum theory of gravity that could describe physical phenomena at the Planck scale (~ 
10 19 GeV/c 2 ), relativists have recently started to explore Lorentz violating theories (LV) (see [I] and references therein). 
Indeed, even though Lorentz invariance (LI) is considered a cornerstone of our knowledge of reality, the challenge 
presented by physics at Planck energy is forcing us to question also our firmest assumptions. In the cosmological context, 
LV theories represent interesting candidates for cosmic acceleration, since in their low-energy limit they generally predict 
a dynamical scalar degree of freedom (DoF) which could provide a source for the late time acceleration, in alternative 
to the cosmological constant. While the standard model of cosmology, based on the laws of General Relativity (GR), 
is to date a very good fit to available data, some outstanding theoretical problems related to the cosmological constant 
have indeed led people to explore alternative theories. To this extent, a wide range of models have been proposed, 
which either introduce a dynamical dark energy (DE) or modify the laws of gravity on large scales (MG) in order to 
achieve self accelerating solutions in the presence of negligible matter. All these alternatives generally result in the 
emergence of new scalar dynamical DoF (see [2, 3]|4] for a comprehensive review), as it is the case with LV theories. 

Interestingly, LV theories typically break LI at all scales, and are therefore constrainable with many different 
measurements and data sets over a vast range of energies. Constraints and measurements on the parameters of a 
general realistic effective held theory for Lorentz violation [5], usually referred to as the Standard Model Extension Eld, 
support LI with an exquisite accuracy. Furthermore, LI has been tested to high accuracy on solar system scales, and 
stringent bounds have been placed on the Post Newtonian parameters (PPN), in particular on those corresponding to 
the preferred frame effects, since such effects are typical of LV theories [Hj - Phenomena on astrophysical scales, and 
in particular tests of gravity in the strong regime, such as those of binary pulsars PEI], provide further bounds on 
LV [S]. On the contrary, the exploration of cosmological bounds on LV theories is still in its infancy 

In the present work, we focus on the class of LV theories known as Hofava gravity [Hill] which modifies the 
gravitational action by adding higher order spatial derivatives without adding higher order time derivatives, thus 
modifying the graviton propagator and achieving a power-counting renormalizability. This is possible if one considers 
that space and time scale differently. Such a prescription is implemented through a breaking of full diffeomorphism 
invariance, which leads to LV at all scales. The resulting theory propagates a new dynamical scalar DoF, i.e. the 
spin-0 graviton. As a candidate for quantum gravity, Hofava theory is expected to be renormalizable and also unitary. 
Nevertheless, at the moment there is no evidence for renormalizability beyond the power-counting arguments. 

Hofava gravity shows a rich phenomenology on cosmological scales, e.g. the higher curvature terms in the action lead 
to a matter bouncing cosmology [1811191 ; it also shows different mechanisms by which it is possible to explain the nearly 
scale invariant spectrum of cosmological perturbations without introducing an inflationary phase [23 Ell 1221 23 E3 , 
finally, cosmological perturbations at late time have been investigated in refs. HU E3 El E3 E3 E3 EDI • 
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In this paper, we perform a thorough analysis of the cosmology in Horava gravity by mapping the theory into the 
framework of Effective Field Theory (EFT) of cosmic acceleration developed in refs. [32 EM EH EH E3 E3 EH] > on the 
line of the EFT of inflation and quintessence [38.1, 39, 40( . The basic idea of this framework is to construct an effective 
action with all the operators which are of relevance to study linear cosmological perturbations around a Friedmann- 
Lemaitre-Robertson-Walker (FLRW) background and are invariant under time-dependent spatial diffeomorphisms. 
Indeed an expanding FLRW background breaks time-dependent diffeomorphism, allowing all these operators to enter 
the action and, furthermore, to be multiplied by a free function of time [38l [35] . The resulting action encompasses 
most models of single scalar field DE and MG which have a well defined Jordan frame. In refs. [32 021 031 04] , the 
EFT framework has been implemented in the public Einstein-Boltzmann solver CAMB [45] 06] , and the associated 
Monte Carlo Markov Chain code CosmoMC 07]. The resulting patches, dubbed EFTCAMB/EFTCosmoMC, are now 
publicly available at http://wwwhome.lorentz.leidenuniv.nl/~hu/codes/ and represent a powerful package which 
allows to explore cosmological constraints both in a model independent and model specific way HU- The original action 
considered in ref. [32], and implemented in the public version of EFTCAMB, contains all Horndeski and some of the 
extensions like GLPV [37] 03] [49] , but does not have all the operators necessary to study Horava gravity. The inclusion 
of Horava gravity in the context of EFT of DE/MG has been recently considered and investigated in refs. 08] I5U] . In 
this paper, we consider the most general action for Horava gravity with all the operators with up to sixth order spatial 
derivatives, which is the minimal prescription to achieve power counting renormalizability. We focus on the part of this 
action that contributes to linear order in perturbations m • For this action we work out a complete mapping to the 
EFT framework deriving also the generalization of the original EFT action used in refs. [32 [32. When we compare the 
predictions of the theory to the observations, we consider only the low-energy operators of Horava gravity, since those 
are the relevant ones to describe the large scale cosmology associated to the observables that we employ. We work out 
the contribution of these operators to the equations of motion for linear scalar and tensor perturbations, implementing 
them in an updated version of EFTCAMB that will be publicly released in the near future. 

The structure of the paper is the following. In section [2j we set up the theoretical background of the paper. In 
particular, in section | 2 . 1 [ we introduce Horava gravity and its main features, while in section | 2 . 2 [ we summarize the 


EFT framework and its implementation in EFTCAMB/EFTCosmoMC. In section 2.3 we work out the mapping of 
Horava gravity in the EFT language focusing on the low-energy part of the action and leaving the mapping of the 
high-energy part of the action to appendix [A] Finally, in section [2~!j we discuss the requirements that EFTCAMB 
enforces on the scalar and tensor DoFs to prevent instabilities in the theory. In section [3] we study the cosmology of 
Horava gravity, discussing in detail how the model is implemented in EFTCAMB and what are the general effects of the 
modifications on the background and the perturbations . Finally, in section[4]we explore observational constraints from 
several combinations of cosmological data sets. To this extent we consider two cases: the low-energy Horava gravity 
action which is characterized by three constant parameters; a subcase of the latter, that evades PPN constraints and 
is characterized by two parameters. We draw our conclusions in section [5] discussing the main results. 


2 Theory 


In this section we set up the theoretical basis for our analysis. In section [2J| we introduce the main aspects of Horava 
gravity, which is the theory we want to investigate and constrain by using the EFT approach. In section 2.2 we review 


the EFT framework, discussing its implementation in EFTCAMB, which is the Einstein-Boltzmann solver we use to 
perform a thorough investigation of the cosmology of the theory. In section |2.3[ we work out the mapping of the 
low-energy Horava gravity action in terms of the EFT functions. The mapping of the high-energy part of the action 
is discussed in appendix [Aj Finally, in section [274] we present the full set of equations evolved by EFTCAMB and the 
conditions that we impose on the tensor and scalar DoFs to ensure that the theory we are considering is viable. 


2.1 Horava Gravity 

Horava gravity has been recently proposed as a candidate for an ultraviolet completion of GR ushizi- The basic idea is 
to modify the graviton propagator by adding to the action higher-order spatial derivatives without adding higher-order 
time derivatives, in order to avoid the presence of Ostrogradski instabilities [52] • The theory is constructed in such a 
way to be compatible with a different scaling of space and time, i.e. 

[dt\ = [ky z , [dx\ = [fc] _1 , (1) 

where z is a positive integer and k is the momentum. In order to accommodate such a different scaling between space 
and time, the action of Horava gravity cannot still be invariant under the full set of diffeomorphisms as in GR, but it 
can be invariant under the more restricted foliation-preserving diffeomorphisms 

x l —t x l (t, x l ) . 
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t —y t (t), 


( 2 ) 




Therefore, within this approach, space and time are naturally treated on different footing leading to Lorentz violations 
at all scales. The emergence of LV is reflected in modified dispersion relations for the propagating DoFs. From 
a practical point of view, the different behavior of space and time is achieved by picking a preferred foliation of 
spacetime, geometrically described within the Arnowitt-Deser-Misner (ADM) formalism. 

It has been shown that the theory is power-counting renormalizable if and only if z > d, where d indicates the number 
of spatial dimensions, which means that the action has to contain operators with at least 2d spatial derivatives f55l 154] . 
Hence, in a four-dimensional spacetime, d = 3, power-counting renormalizability arguments request at least sixth-order 
spatial derivatives in the action. 

Considering the above arguments, the action of Horava gravity can be written as follows m 

Sh — tv — 77 — [ d^Xy/—g (KijK l i — A K 2 — 2 £A + £7 Z + gaidd + L 4 + Lq) + S m \g^ v , x*]> (3) 

i07TGr# J 

where g is the determinant of the metric g^, 1Z is the Ricci scalar of the three-dimensional space-like hypersurfaces, 
I\ij is the extrinsic curvature, and K is its trace. {A, £, 77 } are dimensionless running coupling constants, A is the 
“bare” cosmological constant, a* = <9,;lnA r where as usual N is the lapse function of the ADA! metric. L 4 and Lq denote 
the Lagrangians associated to the higher-order operators, that contain, respectively, fourth and sixth-order spatial 
derivatives (see appendix |A| for the explicit expressions of their parts that contribute to linear order perturbations). 
These Lagrangians constitute the high-energy (HE) part of the action while the operators preceding them represent 
the low-energy (LE) limit of the theory and are the ones of relevance on large scale. S m is the matter action for all 
matter fields, \i- Finally, Gh is the coupling constant which can be expressed as 


G h = £G 


(4) 


where G is the “bare” gravitational constant. As demonstrated in ref. [55, the solution of the static point-like mass 
in the Newtonian limit gives the relationship between the “bare” gravitational constant ( G ) and the Newtonian one 
( Gn ), he. 


G = G n 



(5) 


Then, the coupling in front of the action reads 


1 

I 67 tGh 


W-vY 


( 6 ) 


where mg = 1/87 tGjv is the Planck mass defined locally. 

Notice that the action of GR is recovered when A = 1, £ = 1 and 77 = 0, and the higher order operators in L 4 and 
Lq are not considered. 

The symmetry of the theory allows for a very large number of operators ~ 0(1O 2 ) in L 4 and Lq. In order to 
limit the huge proliferation of couplings in the full theory, in the first proposal Horava imposed some restrictions, i.e. 
projectability and detailed balance (for the details see refs. [551 1511 EH ESI EH 6D : , 1T|). In the following we will not 
impose any of these limitations to the action ([ 3 ]) and we will consider for L 4 and Lq all the operators which contribute 
to the dynamics of linear perturbations m- 


2.2 Effective Field Theory Framework 

In the effective field theory approach to DE/MG [3T1182] . an action is built in the Jordan frame and unitary gauge by 
considering the operators which are invariant under time-dependent spatial diffeomorphisms. The additional scalar DoF 
representing DE/MG is eaten by the metric via a foliation of space-time into space-like hypersurfaces which correspond 
to a uniform scalar field. At quadratic order, which is sufficient to study the dynamics of linear perturbations, the 
action reads 


Seft = 


I d*xV^g{ ^ [1 + f m R + Hr) c(r) a 2 Sg 00 + ^ (aHg 00 ) 2 ^ a 2 Sg 00 SK* 


2 

2 no a 


HM ( SK^) 2 - dK\ 5K\ + m 2 (r) (<T + n^) d^a 2 g™)d„(a 2 g™) 

} ^ HSg 00 61Z + ... 1 + S m [g^ v , Xm], 


(7) 
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where R is the four-dimensional Ricci scalar, Sg 00 , SKSK 1 '^ and 61Z are respectively the perturbations of the 
upper time-time component of the metric, the extrinsic curvature and its trace and the three dimensional spa¬ 
tial Ricci scalar. Finally, S m is the matter action. Since the choice of the unitary gauge breaks time diffeomor- 
phism invariance, each operator in the action can be multiplied by a time-dependent coefficient; in our convention, 
{fi, A, c, Ml Ml Ml Ml Ml ml , M 2 } are unknown functions of the conformal time, r, and we will refer to them as 
EFT functions. In particular, {ST, c. A} are the only functions contributing both to the dynamics of the background 
and of the perturbations, while the others play a role only at level of perturbations. Let us notice that the above action 
includes explicitly all the operators that in ref. [32] have been considered to be relevant for linear cosmological pertur¬ 
bations since they can be easily related to some well known DE/MG models such as f(R), quintessence, Horndeski, or 
because they have been already studied in the EFT of inflation [35J ;40} [31]. For such operators the corresponding field 
equations have been worked out mm- However, additional second order operators can also be considered, such as 
(61Z) 2 , STZ’ jdTZ' 1 ? as well as operators with higher-order spatial derivatives acting on them, [5315211551 155] • In particular, 
as we will show in appendix [Aj additional operators are needed to describe Hofava gravity in the EFT framework (see 
also [48]). 

As mentioned in the Introduction, action 0 allows to describe in a unified language all single scalar field dark 
energy and modified gravity models which have a well defined Jordan frame. In unitary gauge the extra scalar DoF is 
hidden inside the metric perturbations, however in order to study the dynamics of linear perturbations and investigate 
the stability of a given model, it is convenient to make it explicit by means of the Stiikelberg technique i.e. performing 
an infinitesimal coordinate transformation such that r — > t + 7 r, where the new field 7 r is the Stiikelberg field which 
describes the extra propagating DoF. Correspondingly, all the functions of time in action (I7| are expanded in Taylor- 
series and the operators transform accordingly to the tensor transformation laws [53152. Varying the action with 
respect to the 7 r-field one obtains a dynamical perturbative equation for the extra DoF which allows to control directly 
the stability of the theory, as discussed at length in ref. [13- 

In refs. mm the effective field theory framework has been implemented into CAMB/CosmoMC [45, 46j '47] cre¬ 
ating the EFTCAMB/EFTCosmoMC patches which are publicly available at http: //wwwhome. lorentz. leidenuniv. 
nl/~hu/codes/ (see ref. [33] for technical details). EFTCAMB evolves the full equations for linear perturbations 
without relying on any quasi-static (QS) approximation. In addition to the standard matter components (i.e. dark 
matter, radiation and massless neutrinos), massive neutrinos have also been included [33] • EFTCAMB allows to study 
perturbations in a model independent way (usually referred to as pure EFT mode), investigating the cosmological 
implications of the different operators in action |7|. It can also be used to study the exact dynamics for specific models, 
after the mapping of the given model into the EFT language has been worked out (usually referred to as mapping 
mode). In the latter case one can treat the background via a designer approach, i.e. fixing the expansion history and 
reconstructing the specific model in terms of EFT functions; or one can solve the full background equations of the 
chosen theory. We refer to the latter as the full mapping case. Furthermore, the code has a powerful built-in module 
that investigates whether a chosen model is viable, through a set of general conditions of mathematical and physical 
stability. In particular, the physical requirements include the avoidance of ghost and gradient instabilities for both 
the scalar and the tensor DoFs. The stability requirements are translated into viability priors on the parameter space 
when using EFTCosmoMC to interface EFTCAMB with cosmological data, and they can sometimes dominate over 
the constraining power of data [42] . In this paper we will study the case of Hofava gravity, first describing how it can 
be cast into EFTCAMB via a full mapping, then exploring the effects of the stability conditions on its parameter space 
and finally deriving constraints from different combinations of cosmological data sets. 

2.3 Mapping Hofava Gravity into the EFT approach 

In this section we will work out explicitly the mapping of the low-energy (LE) part of action 0 into the EFT formalism 
described in the previous section. This is the part of the action for which we will explore cosmological constraints. We 
show the mapping for the high-energy (HE) part (L 4 and Lq) in the appendix [A] 

We use the following conventions: (-,+,+,+) for the signature of the metric the background is considered 
FLRW with k = 0; dots are derivatives w.r.t. conformal time, r and R = a/a is the Hubble rate; we will use the 
superscript - 0 ' for the background quantities; finally we define a time-like unit vector, n M as 

d t 

Uu = — /! - . with n„n ' 1 = — 1 , ( 8 ) 

M l^d a td 0 l M 

which corresponds to the convention that we use for the normal vector to the uniform-field hypersurfaces in the EFT 
construction of the action 0 [32133. In conformal time and at second order in perturbations, one has 
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(9) 


^ = ^(l+^°° + ^W 0 ) 2 ), 

= g°» (l + l -a 2 5g m + ^(a 2 ^ 00 ) 2 ) . (10) 

In the following, these relations will be often employed. 

Let us first recall the low-energy action, which can be rewritten as: 

2 p 

S H ,LE = (2 ™° ) J d A xV^g (KijK** - A K 2 + £11- 2£A + ga x a l ) 

2 p 

= J n ° —r / d+xV^g i£R + (1 - OK^Kij + (£- \)K 2 - 2£A + rjata') + boundary terms , (11) 

\2£ — V) J 


where the second line has been obtained by using the Gauss-Codazzi relation (62j. 

In the following, we show how to rewrite every single term of the above action in the EFT formalism described by the 
action |7|, providing the mapping of the Horava gravity parameters into the EFT functions. 

• (fSj)(^-2A) 

Comparing the above expression with the EFT action ([7]), it is straightforward to deduce that these two terms 
contribute to the following EFT functions 


(1 + 11 ) 


2£ 

m-vV 


A = -2 


( 2 £ - v) 


A. 


( 12 ) 


« (j^-A)K 2 

In order to identify the relation between the EFT functions and the Horava gravity parameters we have to expand 
AT 2 up to second order in perturbations as 

K 2 = 2I< (0) K + ( SK ) 2 - A' (0)2 , (13) 

by using I\ = K <0> + 5K. Comparing the above relation with the action ([t]) , it is straightforward to see that the 
last term gives contribution to A(r) and the second one to Mf (r). The first term can be computed as follows [31] 


J d i x^2K^K = 2 j d 4 xV^gK^ (V^n") = -2 J d i x v ^j\7 M A' (0) 

-AT(°) 


= 2 


d 4 X\f^g- 


1 - ](«+”) - 


(14) 


where we have integrated by parts the second line and we have used eq. (10). The last line will give respectively its 
contribution to A(r), c(r) and M% (r). Then summarizing, the corresponding contributions to the EFT functions 
from the K 2 term are 


a (t) = - 

M 2 4 (t) = 


^o(£ - A) 
(2? - V) 


A '(°) 2 - 2 


m§(f - A) AT(°) 
'2(2£~V) a ' 


A'(°) 


M 2 (r) = 


f x m 2 (^-A) K<- 0) 
u (2 £-v) a ’ 


2m ° : tt-A). 


(2£ - v) 


(15) 


. ry-ij 

(2£-77) 

As before, we can expand up to second order in perturbations the above operator and it can be written as 

KijK ij = 2K\f5K ij + + SK tj 5K lj , 


(16) 


where we have used the spatial metric to raise the indices and the extrinsic curvature has been decomposed into 
its background and first order perturbation parts, i.e. A', :j - = K ^ 1 +SI\ i j. Moreover, the first term can be written 
as 


2K- ( - ] SK ij = -2-SK = -2% (aK + 3A), 
lJ a a 2 
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(17) 





















where the term proportional to K can be treated as in eq. (14). Finally, in terms of the EFT functions this 
operator can be written as 


A(t) = —- 


n 2 p( 1-Q 
X - v) 
5(i-0 




5(i-0 


c(r) = - ' \U U-H 2 ) 


X - V)a 2 


Mi(r)= " a °; n \ q 2 ), M| = —2 0 


( 1-0 


2o 2 X — 7 7) 


X-»?) ' 


(18) 


J 2 ^ r l^ a 

Let us first write explicitly cii in terms of perturbations up to second order 

<9 ,:N i ^-t„2„oo 


N 


2 ), 


(19) 


where in the last equality we have used a 2 g 00 = — 1 + a 2 5g 00 and then we have expanded in Taylor series. Then 
we get 


m n 


X - v) 


rjdia = 


4X - ? ?) a 




( 20 ) 


where g l i is the background value of the spatial metric. In the EFT language the above expression corresponds 
to 


ml = 


moV 


4(2£ — v)' 


( 21 ) 


Summarizing, we can map the low-energy action © of Horava gravity in the EFT language at the basis of EFTCAMB 
as follows: 


(1 + 0 ) = 


c(r) = — 


2C 


X-X 


TOn 


: 2 X-^) 


(1 + 2£ — 3A) 




Mo = - 


2mn 


x - 0 


-£A - (1 - 3A + 20 


( 1 - 0 . 


•H 2 H 
2a 2 a 2 


M, = -2- m ° 


X-»?) 

TOgjy 

4X-0’ 


(XA), 


M 2 4 (r) = 


2a 2 X-0 


(1 + 2f - 3A) - X) 


M-f = M 2 = 0, 


( 22 ) 


where we have explicitly written the value of the extrinsic curvature and its trace on a flat FLRW background [|] The 
mapping of the high-energy part of the action can be found in appendix |Aj 


2.4 Degrees of freedom: dynamics and stability 


After the full diffeomorphism invariance is restored by means of the Stiickelberg mechanism, at the level of perturbations 
we have a dynamical equation for the scalar DoF represented by the Stiickelberg field n. In the case of the low-energy 
limit of Horava gravity that we are considering, this equation reads 


r]TT + 2ri'H'k + 


(3A - 2f - 1)X 2 — H)+ riCH 2 + U) ir + k 2 £(\- 1)tt + £(A- l)kZ + 


(Xi)X-q) 

2 k 


v a 2 (pi +pi) 


~Vi 


= 0, 
(23) 


1 For the low-energy action it is possible to obtain part of the mapping by following the method in ref. m- However, one has to consider 

that our formalism and notation differ from the one in ref. m because we are using conformal time, a different signature for the normal 

unit vector, a different notation for the EFT functions and one more operator is included in our low-energy action: a^a^. 
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where Z is the standard CAMB variable mmm Pi, pi are the background density and pressure of matter com¬ 
ponents, and Vi is the velocity perturbation of matter components. The above equation is coupled with the following 
perturbative field equations: 

• time-time (t) held equation 

2 n [fc 2 7r(?7 - 3A + + 1) + (1 - 3A )kZ] + 2 k 2 (2£fj + pn) + a 2 ^^-5 Prn = 0 , (24) 

m 0 


• space-space (s) held equation 

-AH \k 2 ( 3A - 2£ - 1)tt + (3A - l)kZ] + (1 - 3A )h + 4fc 2 ^ + 2fc 2 (-3A + 2£ + 1)tt + 3a 2 ^ ~ ^ 6P m = 0, (25) 

TOq 

where h, fj are the usual scalar perturbations of the metric in synchronous gauge (notice that we have added a bar to 
the standard metric perturbation in order to do not confuse it with the Horava gravity parameter, rj). EFTCAMB 
evolves the above set of coupled differential equations along with the usual matter perturbation equations and the 


initial conditions are set following ref. m- Let us notice that by using the mapping \22\ worked out in the previous 
section, it is straightforward to deduce the above equations following the general prescription in ref. [5T] . 

We shall now determine the dispersion relation of the scalar DoF, computing the determinant of the matrix of the 


coupled system eqs. (23)- (25). Since the number counting of dynamical DoFs will not be changed by neglecting the 


couplings with standard matter species, for simplicity, for the purpose of this calculation we neglect them. After taking 
the Fourier transform d T —> —ito, we can rewrite the system (23)- (25) in the following matrix form: 


0^777r 

7Vh 

7tT7? 

^STT 

7sh 

7s fj 

"7t7T 

7th 

Itfj 



(26) 


where the term j a b with a,b = {ir, h, r)} corresponds to the coefficient of b in equation a and they can be easily deduced 
from the above equations. Finally we set the determinant to zero and get 


k 4 to(uj + iH) 


to 2 + i 2 Hto — 


(A — 1)£(2£ — rj) 2 Z 


»?(3A-1) 

which can be written in a compact form as 


AT - - \{H- - 3A + 2£ + 1) + (6A - - 2)H‘ 


krui 


+ i^j [to 2 + iato — fc 2 c 2 + 0 \ = 0 . 


= 0 , (27) 


(28) 


From the above equation we deduce that only one extra dynamical DoF exists, which corresponds to the scalar graviton 
(tt field in EFT language), as expected. Furthermore, one can identify the terms in the squared bracket as follows: a 
is a friction term, /? is the dispersion coefficient and c 2 can be identified with the canonical speed of sound defined in 
vacuum, when no friction or dispersive terms are present. Let us notice that both the friction and dispersive terms 


are related to the nature of the dark energy component through the dependence of the Hubble rate on the latter (381. 


but here we include also friction and 


The procedure to compute the dispersion relation (27) follows the one in ref. 
dispersive terms. 

In order to ensure that a given theory is viable, we enforce a set of physical and mathematical viability conditions. 
The mathematical conditions prevent exponential instabilities from showing up in the solution of the 7r-field equation, 
and the physical ones correspond to the absence of ghosts and gradient instabilities for both scalar and tensor modes. 
In particular, in our analysis of Horava gravity, for the scalar DoF they correspond to 


2toq77(1 — 3A )k 2 


(t? - 20 (2(3A - 1 )U 2 + t?(A - 1 )k 2 ) 


> 0 , 


Q2C-0(A- 1) 

7y(3A - 1) 


> 0 , 


(29) 


where the first condition corresponds to a positive kinetic term and it has been obtained from the action by integrating 
out all the non dynamical fields, while the second one ensures that the speed of sound is positive. Let us note that the 
ghost condition reduces to the one in the Minkowski background by setting the limit a —> 1. 

Additional conditions to be imposed comes from the equation for the propagation of tensor modes hij, 


A T {r)hij + B t (t )hij + D T (T)k 2 h i:j + E Tij = 0. 


(30) 















where STij generally contains the matter contributions coming from the neutrino and photon components and, for 
Horava gravity, the remaining coefficients read: 


At = 


2 


Dt 


2 £ 

2 ? - r, ’ 


Bt = 


E' 


Tij 


4 n 

2 Z-r,’ 

- —XT 

~ ml lJ ' 


(31) 

(32) 


The viability conditions require At > 0 and Dt > 0 to prevent respectively a tensorial ghost and gradient instabili¬ 
ties [ 33 ] . 

It is easy to show that the above conditions translate into the following constraints on the parameters of Horava 
gravity: 

0 < 77 < 2 £, A>1 or A < -, (33) 

o 

which are compatible with the viable regions identified around a Minkowski background m- In the following we will 
not explore the A < 1/3 branch since along it the cosmological gravitational constant on the FLRW background becomes 
unacceptably negative [B2'i lS 3( and the branch does not have a continuous limit to GR. The conditions that we have 
discussed are naturally handled by EFTCAMB/EFTCosmoMC in the form of viability priors that are automatically 
enforced when the parameter space is being sampled. 

Besides the above theoretical viability conditions, there are observational constraints on the Horava gravity param¬ 
eters coming from existing data. In particular: 


Big Bang Nucleosynthesis (BBN) constraints [TT] , which set an upper bound on \G cosmo /GN — 1| < 0.38 (99.7% 
C.L.fl where G cosmo is the cosmological gravitational constant as defined in section 


3.1 


Solar system constraints, where the parametrized post Newtonian parameters (PPN) are bounded to 

a, < 3.0 ■ 10 -4 (99.7%C.L.), a 2 < 7.0 ■ 10 “ 7 (99.7%C.L.). (34) 


where a\ and oli are two of the parameters appearing in the PPN expansion of the metric around Minkowski 
spacetime, more precisely those associated with the preferred frame effects mm ■ Here we consider only these 
two parameters since they are the only ones of relevance for constraining LV. It has been shown in refs. [6611671 58] , 
that the PPN parameters for the low-energy action of Horava gravity, read 


= 4(2£ -r,- 2 ), 

(t? - 2£ + 2)( 77 ( 2 A - 1) + A(3 - 4£) + 2£ — 1) 
(A-l)(„-20 


(35) 


It is easy to show that combining the above relations, the above mentioned PPN bounds result in a direct 
constraint on A that reads: 


log 10 (A - 1) < -4.1 (99.7%C.L.), (36) 

while the bound on au provides a degenerate constraint on the other two parameters {■f, 77 }. 

• Cherenkov constraints from the observation of high-energy cosmic rays [J>Sj are usually imposed as a lower bound 
on the propagation speed of the scalar DoF and the propagation speed of tensor modes. In the case of LV theories 
we will refer the reader to Refs. EH [70], for further details. However, since these bounds have not been worked 
out specifically for Horava gravity we decided not to impose them a priori. 

For the present analysis we consider two specific cases of Horava gravity: 

1 . Horava 3, hereafter H3, where we vary all three parameters {A, 7 y,£} appearing in the low-energy Horava gravity 
action; 

2 The original bound in ref. m is reported at 68% C.L. and we convert it to 99.7% C.L. by assuming a Gaussian posterior distribution 
of Gcosmo /G n — 1. 

3 The original bounds in ref. EH (and references therein) are reported at 90% C.L. and we convert it to 99.7% C.L. by assuming a 
Gaussian posterior distribution of the relevant parameters. 
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2. Hofava 2, hereafter H2, where we choose the theory parameters in order to evade the PPN constraints (35) by 
setting exactly a\ = a .2 = 0. This implies: 

V = 2£ - 2, (37) 


so that the number of free parameters reduce to two, {A, 77 }. This case has the quality of systematically evading 
solar system PPN constraints, meaning that it is not possible to build a local experiment, with arbitrary precision, 
to distinguish it from GR. Therefore it can only be constrained with cosmological observations. 


For both cases we impose the physical and mathematical viability conditions in the form of viability priors as 
discussed in ref. |42| . The portion of the parameter space excluded by the viability priors can be seen as a dark grey 
contour in figure [6] for the H3 case and in figure [7] for the H2 case. For both cases we also derive the bounds on 
Gcosmo/GN — 1 and for the H3 case we provide cosmological bounds on the PPN parameters. These results are shown 
and discussed in detail in section [4j 


3 Hofava Cosmology 

In this section we highlight the cosmological implications of the low-energy Hofava gravity cases, H2 and H3, previously 
introduced. In section |3.1| we discuss the changes that Hofava gravity induces at the level of the cosmological back¬ 
ground, while in section |3.2| we elaborate on the effects that are displayed by the theory at the level of perturbations 
by means of two examples. 


3.1 Background 


The first step towards testing a theory against cosmological observations, is to investigate the behaviour of its cosmolog¬ 
ical background. In this section, we discuss the background evolution equation for Hofava gravity, its implementation 
in EFTCAMB, and review the definitions that we adopt for the cosmological parameters. 

The Hofava gravity field equations for a flat FLRW background read: 


3A - 1 


2 

3A - 1 


H = 


87 rG N ( 2 £ - rj) 


A 


6 


[ J2pi 


1 


H + -W 




= -- 7 —a" + 47 tGn 


( 2 £ - v) 




-a ? Vi 


(38) 

(39) 


where pi and pi are respectively the density and the pressure of the matter fluid components, i.e. baryons and dark 
matter (m), radiation and massless neutrino (r) and massive neutrinos ( v ). In this work we consider that all massive 
neutrino species have the same mass and we set the sum of their masses to be 0.06 eV. In addition to the Friedmann 
equations, we have the standard continuity equations for matter and radiation: 


Pi + 37f(l + Wi)pi = 0, 


(40) 


while for massive neutrinos we refer the reader to ref. l43l for a detailed discussion. 


Starting from the Friedmann eq. (38), we can define the cosmological gravitational constant as: 


G c 


(2g ~ V) 
3A — 1 


G n , 


(41) 


differs from Gn, which is obtained with local experiments, as already pointed out in 


where it is clear that G„ 
ref. [511 . This definition allows us to write the Friedmann equation (38) in another way: 


U = 87 rG c 


'Yhi Pi 


1 


2£ A 


87 tG n 2£ - rj 3 


(42) 


From this equation it is straightforward to see that in general, once the theory parameters have been properly set, the 
modification that Hofava gravity induces at the level of the background is a global rescaling of H Ca¬ 
in order to properly identify the parameters that we should fit to data, we have to pay special attention to the 
working definition of all the relevant quantities. In particular in the definition of the relative density abundance. For 
the matter fields, we define f]j(a) in terms of the locally measured gravitational constant, Gn, and the present time 
Hubble parameter, Hq. We then derive the abundance of the effective dark energy, describing the modifications to the 
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Figure 1: The figure shows the evolution of the densities parameters for baryons and dark matter (m, dashed line), 
radiation, neutrino and massive neutrinos (r, dot dashed line) and dark energy (DE, solid line). In the left panel we 
compare the density parameters of the H3 case (green lines) with the ones in the ACDM model (red lines). In the 
right panel the comparison is between the H2 case (blue lines) and ACDM. The yellow area highlights the radiation 
dominated era. For this figure the standard cosmological parameters are chosen to be 12° h 2 = 0.0226, 12° h 2 = 0.112, 
12° h 2 = 0.00064 and H 0 = 70Km/s/Mpc. In the H3 case the Horava gravity parameters are A = 1.4, £ = 0.9, = 1.0 

while in the H2 case they are fixed to A = 1.4, r] = 1.0. 


Friedmann equations, by means of the flatness condition, i.e. ]TT ^*( a ) + ^de(o) = 1- To this extent, we rewrite the 
Friedmann eq. (38) as 

,2 o„r< Si Pi 2 i 2£ K * . f - 3A — 1 


H = 8t tG n - 


so that it is straightforward to identify 


2£-p3 


a + 1 — 


2£ — V 


H 


(43) 


12, (a) = 

n , , 2£ A a 2 1 3A - 1 

D_d_e(o) = —-— —^ + 1 - —-■ 

2£ — r] 3 2£ - rj 

At present time (ao = 1), we can immediately see that = 1 — 12° with: 


0 _ 2£ A 

DE 2^ — 77 3 H 2 


+ 1 - 


3A - 1 
2£-p' 


This allows us to rewrite the Friedmann eq. (38 1 in terms of the parameters that we are going to sample as: 


7i = 


(2£ - 7?) 
3A — 1 


a z H 2 


0 ° 

a L m 

a 3 


$2^ 

a 4 


^DE ~ 1 


3A - 1 
2 


(44) 


(45) 


(46) 


This is the background equation that EFTCAMB evolves, along with its time derivatives. For details about how the 
code treats see ref. [43] . Finally, one can use eq. (451, to substitute the “bare” cosmological constant with fl° DE , 


therefore in the following we use the latter as one of the Horava parameters that we fit to data instead of A. 

We shall now specialize to some choices of the Horava parameters, and derive the corresponding expansion history 
in order to visualize and discuss the effects of Horava gravity, in particular for the H3 and H2 cases, on background 
cosmology. We choose the background values of the cosmological parameters to be 12° h 2 = 0.0226 for baryons, 
f2° h 2 = 0.112 for cold dark matter, 12° h 2 = 0.00064 for massive neutrinos and Hq = 70Km/s/Mpc, accordingly to the 
default CAMB parameters. Additionally, the parameters of the H3 case are chosen to be: A = 1.4, £ = 0.9, 77 = 1.0; 
while in the H2 case we set A = 1.4 and 77 = 1.0. While the general trend of the modifications does not depend on the 
magnitude of the theory parameters, the above values are selected in order to enhance the effects and clearly display 
the changes with respect to the standard cosmological model, ACDM. Thus they have to be considered as illustrative 
examples because the values involved are significantly bigger than the observational bounds that we will derive in 
section [3] However, in both cases the choices of parameters respect the viability criteria discussed in section [24] 
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In figure[T]we can see the behaviour of the relative densities for matter (dark matter and baryons), radiation (photons 
and relativistic neutrinos), and effective dark energy, as defined in eqs. (44). One can notice that at early times the 
matter species display density values that are generally bigger than one, on the contrary the dark energy component 
assumes negative values. This can be explained as follows. The matter components are well behaved, with positively 
defined densities with a time evolution that is exactly the standard one (eq. ( |40[ ) ), as expected when working in Jordan 
frame. However, the expansion history changes as it is rescaled by a constant (eq. ©), altering the time behaviour of 
the relative abundances. The effective dark energy balances this effect in order to respect the flatness condition. We 
argue that in this specific case the interpretation of the modification of gravity in terms of a fluid-like component is not 
well justified/posed, representing instead a genuine geometrical modification of the gravitational sector. This kind of 
behaviour for the effective dark energy component is commonly encountered in dynamical analysis studies of modified 
gravity models where the flatness condition is used as a constraint equation EDGE]- From figure [I] we can also notice 
that Horava gravity does not affect the time of radiation-matter equality as the continuity equations for these species 
are not changed, as it is clearly highlighted by the yellow region in the figure. Indeed and f l r for all the models 
cross at the same value of the scale factor. On the other hand, the time of equality between matter and dark energy is 
slightly modified depending on the model parameters. Finally, let us notice that, once the parameters of the theory are 
chosen to be compatible with the observational constraints, all these effects that we have discussed are quite mitigated 
and become hardly noticeable by eye in the plots. Indeed, values of the parameters consistent with the bounds that 
we derive in section [4] would induce a less negative DE density at earlier times. 


3.2 Perturbations 

In this section, we proceed to study the dynamics of cosmological perturbations. Once we have worked out the 


background equations of Horava gravity (46), as well as the mapping of this theory into the EFT language (22), we 


have all the ingredients required by EFTCAMB to perform an accurate analysis of the perturbations. For technical 
details on the actual implementation, as well as the full set of perturbative equations that are evolved by EFTCAMB, 
we refer the reader to ref. [33j. 

As we will see, the behaviour of perturbations in Horava gravity displays an interesting and rich phenomenology, 
allowing to investigate the theory and to constrain its parameters with the available data. In the following, we perform 
an in depth analysis of the dynamics of linear perturbations and the corresponding observables, specializing to a choice 
of parameters for the case H3 and one for the case H2, in order to visualize and quantify the modifications. In all cases, 
we set the values of the cosmological parameters to the one used in the previous section, which are the default CAMB 
parameters, while for the Horava parameters we use: in the H3 case, (£ — 1) = —0.01, (A — 1) = 0.004, p = 0.01; in 
the H2 case, (A — 1) = 0.02, p = 0.05. As it will be clear in the next section, these are noticeably bigger than the 
observational constraints that we will derive, but they facilitate the visualization of the effects on the observables. Let 
us stress that, while the direction and entity of the modifications that will be described in the remaining of this section 
are specific to the choice of parameters, we have found an analogous trend for several choices of parameters that we 
have sampled in the region allowed by the viability priors. 

Let us now focus on the time and scale evolution of cosmological perturbations and the growth of structure. In 
order to discuss the deviations of Horava gravity from ACDM, we study the behaviour of the p(k, a)-function, which 
is defined in Newtonian gauge as |72] 


k 'F = -p(k,a)—^p m A 


2toq 


(47) 


where A m is the comoving matter density contrast and *F is the scalar perturbation describing fluctuations in the time¬ 


time component of the metric. As it is clear from eq. (47), p parametrizes deviations from GR in the Poisson equation. 


In the standard cosmological model, ACDM, this function is constant and p = 1. Let us notice that EFTCAMB does 


never evolve the above quantity (47), but it can easily output /i as a derived quantity. Moreover, we also analyse the 
behaviour of the quantity <F + *F, where $ is the scalar perturbation of the space-space component of the metric in 
Newtonian gauge. This quantity is important as it allows to identify possible modifications in the lensing potential and 
in the low multipole of the cosmic microwave background (CMB) radiation through the Integrated Sachs-Wolfe (ISW) 
effect. Finally, we explore the fluctuations in the total matter distribution defined as 6 m = )T) m p m A m / p m . 

In figure [2] we show the time and scale behaviour of these three quantities. In order to facilitate the visualization of 
the deviations from the ACDM behaviour, we show the logarithmic fractional comparison between these quantities in 
the two Horava gravity cases considered and the ACDM model. 

• H3 case: from the top left panel of figure [2] we can see that p significantly deviates from one at large scales and 
all redshift with fractional differences that are around unity (100%). Small deviations of the order of 10 -4 can 
be also seen at small scales and high redshift. At small scales and low redshift, in the bottom right part of the 
H3 a) panel, one can notice small features due to the fact that the n field oscillates while being coupled to the 
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Figure 2: We show the relative comparison of the modification of the Poisson equation p, the source of gravitational 
lensing $ + U/ (whose derivative sources the ISW effect on the CMB), and <5 m = p m A m / J2 m Pm with their ACDM 

values for the H3 (upper panel) and H2 (lower panel) models. In all panels, the dashed white line represents the physical 
horizon while the solid white line shows where the relative comparison changes sign. For this figure the standard 
cosmological parameters are chosen to be 12° h 2 = 0.0226, h 2 = 0.112, f2° h 2 = 0.00064 and Hq = 70Km/s/Mpc. In 
the H3 case the additional parameters are (£ — 1) = —0.01, (A — 1) = 0.004, r) = 0.01 while in the H2 case they are 


fixed to (A — 1) = 0.02, rj = 0.05. For a detailed explanation of this figure see section 3.2 


other species. From the top central panel of the same figure we can see that gravitational lensing is modified 
as well. On large, super-horizon, scales deviations from the ACDM behaviour are not significant, staying below 
10~ 2 at all the times shown. In general at these scales the lensing is suppressed. On sub-horizon scales in turn 
the enhancement of the lensing potential with respect to the ACDM case becomes relevant. A similar behaviour 
can be seen in the total matter density contrast. Although on super-horizon scales, as well as just below the 
horizon, the density contrast is enhanced compared to the ACDM one, on very small scales it is suppressed. 
Noticeably the oscillations that we see in p do not reflect on $ + 'P and S m , which look rather regular. The 
physical interpretation of this is that even if the additional scalar DoF is introducing fluctuations in the structure 
of the Poisson equation the held is not coupled strongly enough to introduce fast fluctuations in the matter and 
metric fields themselves. 

• H2 case: from the lower left panel of figure [2] we can notice that, in the H2 case, the behaviour of p is rather 
different from the H3 case. In particular on small scales its value returns to the GR one. This is compatible 
with the extra constraint that we have imposed in this case ( |37[ ), making the theory indistinguishable from GR 
on solar system scales. On large scales and high redshift, similarly to the H3 case, deviations from the ACDM 
behaviour are of the order 10 -2 (1%). Panels H2 b) and H2 c) in figure [ 2 ] show that the lensing effects and the 
growth of matter perturbations do not follow the trend of p. Indeed, in the case of lensing, in panel H2 b), around 
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Figure 3: We show the quantity £/v = 'Kn/'Httn that we introduce as an indicator of the goodness of the quasi-static 
approximation for the H3 (left panel) and H2 (right panel) cases. In both panels, the dashed white line represents the 
physical horizon, while the solid white lines highlight the scale dependence of this quantity. For this figure the standard 
cosmological parameters are chosen to be 12° h 2 = 0.0226, 12° h 2 = 0.112, 12° h 2 = 0.00064 and Hq = 70Km/s/Mpc. In 
the H3 case the Horava gravity parameters are (£ — 1) = —0.01, (A — 1) = 0.004, rj = 0.01 while in the H2 case they 


are fixed to (A — 1) = 0.02, rj = 0.05. For a detailed explanation of this figure see section 3.2 


and below the horizon the model displays significant deviations from the ACDM behaviour that are similar to 
the H3 ones. From panel H2 c) we notice that the growth of matter perturbations deviates significantly from the 
ACDM one (around 10 -2 ) at almost all redshifts and scales. Finally, in the same panel it can be noticed that the 
density contrast is enhanced for k < 10 -1 h/Mpc while it is suppressed at very small scales and all redshift. 


After considering the cosmological evolution of metric and matter perturbations we now turn to the study of the 
dynamics of the additional scalar DoF that propagates in Horava gravity. In particular we study the quantity introduced 
in ref. [41] to quantify the deviations from quasi-staticity for the dynamical scalar DoF, 7r, i.e. 


6v = 


T~Lttn 


(48) 


where with the index N we indicate that we are working with the 7r-field in Newtonian gauge. This quantity compares 
the evolutionary time-scale of the additional scalar DoF with the Hubble time scale, thus quantifying how many times 
the 7r-field changes significantly in a Hubble time. Small values of this quantity imply that the 7r field is slowly evolving 
and that time derivatives of the field can be neglected when compared to the value of the scalar field itself. On the 
contrary large values mean that the time derivative of the field is playing a major dynamical role, and hence QS would 
not be a safe assumption. 

The time and scale behaviour of £at can be seen, for the H3 and H2 cases, in figure [3j We can notice that, roughly 
for both cases, the 7r-field is slowly evolving at low redshift (0 < z < 1), on the other hand, at higher redshift we can 
see that its dynamics becomes relevant and deviations from a QS behavior are order 30%. We can also notice that, at 
all scales and times, the evolutionary time scale of 7 tjv is smaller than the Hubble rate. From the same figure we can 
see that this evolutionary rate does not significantly depend on scale. The white lines in figure [3] show some residual 
scale dependence at early times and clearly show that this scale dependence gets weaker at late times. 

Finally, we discuss how the modified dynamics of perturbations in Horava gravity affects the observables that we 
later use to constrain this theory. In figure |4j we compare several power spectra for the H2 and H3 cases in comparison 
to the ACDM model. We identify the following effects on the observables: 

• Differences in the late time Integrated Sachs-Wolfe (ISW) effect. For the two cases that we explore, we find an 
enhancement of the amplitude of the low-£ temperature power spectrum, as it can be seen from the top left panel 
of figure [4] which is related to an increase of the late-time ISW effect m- The latter is sourced by the time 
derivative of + 'F and, as we can see from figure [2j for the two Horava gravity cases the time evolution of this 
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Figure 4: Power spectra of different cosmological observables in the ACDM, H2 and H3 cases. Upper panel: CMB 
temperature-temperature power spectrum at large (left) and small (right) angular scales. Central panel: lensing 
potential and CMB temperature cross correlation power spectrum (left), lensing potential auto correlation power 
spectrum (right). Lower panel: matter power spectrum (left) and B-mode polarization power spectrum (right). In this 
last panel the solid line corresponds to the scalar induced B-mode signal while the dashed one shows the tensor induced 
component. For this figure the standard cosmological parameters are chosen to be f l®h 2 = 0.0226, 52° h 2 = 0.112, 
52° h 2 = 0.00064 and H 0 = 70 Km/s/Mpc. The Horava gravity parameters in H3 case are chosen to be: (£— 1) = —0.01, 
(A — 1) = 0.004, r) = 0.01; in the H2 case they are: (A — 1) = 0.02, rj = 0.05. 


quantity is modified. This change also affects the CMB temperature-lensing cross correlation (central left panel), 
as discussed below. 

• Differences in the gravitational lensing. As we already discussed, in the specific cases that we explore, gravitational 
lensing results to be enhanced as we can see in the central panel of figure [2j This reflects on the CMB lensing 
power spectrum as shown in the central right panel of figure [4j where we can notice that fluctuations of this 
observable are enhanced for both H3 and H2 cases with respect to the ACDM model. This modification also has 
an effect on the high multipole of the lensed CMB temperature power spectrum as highlighted in the top-right 
panel of figure [4] At first glance we can see that, compared to the ACDM model, the profile of the high-5 peaks is 
less sharper in the H3 and H2 cases because of the lensing enhancement. We can also notice that there is a slight 
asymmetry between peaks and troughs due to a combined effect of the lensing modification with the modified 
Hubble rate discussed in section [3J| thus leading to a small change in the angular scale of the CMB peaks. From 
the central left panel, we can see that the CMB temperature-lensing cross correlation spectrum is influenced by 
both the ISW and lensing modifications. In particular, this spectrum results to be enhanced at low-5 because of 
the lensing and ISW enhancements but it is suppressed for 50 < £ < 100 following the trend of the temperature 
power spectrum. Indeed, we can notice, from the top right panel of figure [4] that at these scales the spectra are 
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Figure 5: The tensor induced component of the B-mode CMB polarization power spectrum in the ACDM, H2 and 
H3 cases. For this figure the standard cosmological parameters are chosen to be 0° h 2 = 0.0226, h 2 = 0.112, 
h 2 = 0.00064, r = 1 and H 0 = 70Km/s/Mpc. The Horava gravity parameters in H3 case are chosen to be: 
(£ — 1) = —0.3, (A — 1) = 4 x 10~ 4 , g = 10 -3 in the H2 case they are: (A — 1) = 1, rj = 0.6. 


suppressed due to the lensing effect as previously mentioned. Finally, the enhancement of the lensing potential 
also affects the component of the CMB B-mode power spectrum that is sourced by the lensing of the E-mode of 
polarization. This situation is highlighted in the lower right panel of figure [3] The solid lines representing this 
component of the B-mode spectrum are enhanced proportionally to the enhancement in the lensing potential. 

• Differences in the growth of matter perturbations and the distribution of the large scale structure. For the two 
cases under analysis (H3 and H2), we observe a slight enhancement of the growth of structure in the total matter 
power spectrum, at intermediate scales, as well as a slight suppression on small scales, as it is clearly depicted in 
the lower left panel of figure [fj and in agreement with our previous analysis of the density contrast, see figure [2] 
The matter power spectrum, for both H2 and H3, follows the ACDM one on large scales (fc < ICC 3 h/Mpc) while 
for 10 -3 < k < 1CT 1 h/Mpc it is slightly enhanced, particularly for the H2 case. At very small scales, both the 
H2 and H3 matter spectra follow the ACDM behaviour. 

• Differences in the propagation of tensor modes. As previously discussed in section |2.4[ the tensor dynamical 
equation is modified in Horava gravity. This change is usually reflected in the tensor induced component of 
the B-modes of CMB polarization m Eg. In particular, in the H3 case, the parameter £ controls directly 
the propagation speed of gravitational waves, while the combination 2£ — is responsible for the strength of 
coupling between tensor modes and matter. Instead, in the H2 case the tensor speed of sound is controlled by rj, 
while there is no effect on the coupling with matter. The choice of parameters we made for figure [4] displays a 
significant effect on the scalar component of the B-mode spectrum as shown in the lower right panel of figure [4] 
as solid lines, but the effect on the tensor component (dashed line) of the B-mode power spectrum for the same 
parameters is much smaller and not visible in the figure. In figure [5] we change the Horava gravity parameters to 
better display the effect of the change in the tensor sector. Therefore only for this figure we choose the Horava 
gravity parameters in H3 case to be (£ — 1) = —0.3, (A — 1) = 4 x 10 —4 , 77 = 10 -3 and in the H2 case they are 
(A — 1) = 1,77 = 0.6. As we can see from that figure, the leading effect is due to the modification of the speed of 
gravitational waves [75]. In the next section we will find that due to a combination of viability requirements and 
data constraints, for the H3 case, £ < 1, therefore the spectrum results to be shifted to the right with respect to 
the ACDM one, since tensor modes propagate sub-luminally. On the other hand, in the H2 case, tensor modes 
propagate super-luminally {p > 0) and the whole spectrum is shifted to the left. Finally, a modification of the 
coupling to matter leaves an observational imprint that is much smaller that the previous one as cosmological 
gravitational waves propagate almost in vacuum. We can also conclude, on the basis of the results we will present 
in the next section, that in the H3 case since the tensor sound speed is less than one, the Cherenkov constraints 
are not always satisfied but only in a very tiny range mm- On the contrary, in the H2 case the tensor sound 
speed is always super-luminal, then the Cherenkov constraints are evaded. 


4 Cosmological constraints 

In this section we derive and discuss the observational constraints on Horava gravity coming from cosmological probes. 
After describing the data sets used, we focus on the H3 and H2 cases described at the end of section [2~T[ 
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1 + 

O C 

b c 

to I- 

-0 01+ U - Ui 

U.U1—0.02 

-0 01+°- Oi 

u - ul -0.02 

-0 01+ U Ui 

U ' U1 -0.02 

l°gio(A - 1) 

[-8,-2] 

< -4.56 

< - 4.68 

< -4.24 

< -4.68 

< -4.31 

logio V 

[-5,-1] 

— 

— 

— 

— 

— 

no 

lL DE 

derived 

0.69t°;° 6 7 

0.69 ±0.02 

Q 7Q+0.04 

u - ' u -0.05 

0.66 ± 0.06 

0.69 ±0.02 

otx 

derived 

< 0.283 

< 0.240 

< 0.263 

< 0.322 

< 0.220 

logio «2 

derived 

< 5.05 

< 4.83 

< 5.04 

9 '7D+2.52 
z *' u -2.41 

< 4.72 

^cosmo /Gn 1 

derived 

< 0.035 

< 0.030 

< 0.033 

< 0.040 

< 0.028 


Table 1: The 99.7% C.L. marginalized posterior bounds on the H3 case parameters and relevant derived parameters. 


4.1 Data sets 

In our analysis we use several geometrical and dynamical probes, combining them progressively. 

The first data set employed, hereafter PLC, consists of the low-f (2 < £ < 50) CMB temperature-temperature 
power spectra from the Planck satellite [7i[73, considering the 9 frequency channels ranging from 30 ~ 353 GHz. 
In addition, we consider the 100, 143, and 217 GHz frequency channels for the high-£ modes (50 < £ < 2500) of the 
CMB temperature spectrum. We also include the WMAP \ow-£ polarization spectra (2 < £ < 32) 7S] to break the 
degeneracy between the re-ionization optical depth and the amplitude of CMB temperature anisotropy. 

The second data set considered, denoted as BG, is a combination of background measurements that helps to break 
the degeneracies between background parameters and the ones describing the behaviour of perturbations. We use data 
from HST [79] which measures the local Hubble constant from optical and infrared observations of more than 600 
Cepheid variables. In addition, we consider the “Joint Light-curve Analysis” (JLA) Supernovae sample as analysed in 
ref. [80] which is constructed from the SNLS, SDSS and HST SNe data, together with several low redshift SNe. We 
also employ baryon acoustic oscillations measurements taken from: the SDSS Main Galaxy Sample at z e g = 0.15 [5Tj : 
the BOSS DR11 “LOWZ” sample at z eS = 0.32 [52]; the BOSS DR11 CMASS at z eff = 0.57 of [52]; and the 6dFGS 
survey at z e g = 0.106 [53] . 

The third data set that we use, consists of the Planck 2013 full-sky lensing potential power spectrum obtained by 
using the 100, 143, and 217 GHz frequency bands employed by the Planck collaboration to detect the CMB lensing 
signal with a significance greater than 25<r [84]. We refer to this data set as the lensing one, hereafter CMBL. 

Finally, we use the measurements of the galaxy power spectrum by the WiggleZ Dark Energy Survey [55] in order 
to exploit the constraining power of large-scale structure data. The WiggleZ data set consists of the galaxy power 
spectrum measured from spectroscopic redshifts of 170,352 blue emission line galaxies over a volume of 1 Gpc 3 [HIST]. 
The covariance matrices are taken to be the ones given in ref. m and are computed using the method described in 
ref. [90] . It has been shown in refs. m S3 that linear theory predictions are a good fit to the data regardless of 
non-linear corrections up to a scale of k ~ 0.2h/Mpc. Since changes in the growth induced by modifications of gravity 
can slightly alter this scale, in this work, we use the WiggleZ galaxy power spectrum with a more conservative cut of 
fc max = O.lh/Mpc. We marginalize over a scale independent linear galaxy bias for each of the four redshift bins, as in 
ref. [87] . Let us notice that, in general, the linear galaxy bias in a modified gravity scenario is not scale independent, 
as shown in refs. 1851189] . However, as we will show in the next section, the cosmological constraints are mainly driven 
by CMB and background observables towards the ACDM limit. Therefore, we do not expect that a scale dependent 
bias will dominate the results in Horava gravity when considering WiggleZ data. 


4.2 H3 case: results 


The first case we compare to cosmological observations is the low-energy limit of Horava gravity, H3, for which the 
parameters of the theory are {£, A, 77 }. 

We sample A and £ shifting them by one so that the GR limit of the new parameters, i.e. A — 1 and £ — 1, corresponds 
to a zero value. In addition, we use a logarithmic sampler for the parameters i] and A — 1 as they are constrained to 
be positive by physical viability, as discussed at length in section [2~T| Instead, £ is allowed to change sign, hence we 
sample £ — 1 linearly. 


When combining the viability priors discussed in section 2.4 with cosmological data, we notice that the requirement 
of physical viability has a strong effect on the posterior distribution of the parameter £. This is shown in panel (a) of 
figure [ 6 ] Even though £ is not constrained to be above or below 1 a priori, the magnitude of A, draws the posterior of 
the model into a region where only values of £ < 1 are viable. Let us notice that there is a very small viable region 
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Figure 6 : Results of our analysis of the H3 case. Upper panel: The marginalized joint likelihood for combinations of the 
parameters of low-energy Horava gravity. The darker and lighter shades correspond respectively to the 68 % C.L. and 
the 95% C.L.. Lower panel: The marginalized likelihood of the parameters of low-energy Horava gravity. In both panels 
different colours correspond to different data set combinations as shown in legend. The dark grey shade corresponds 
to the marginalized region of parameter space excluded by viability priors. 


in the parameter space above one (1 < £ < 10 6 ), which due to the scale adopted in the plot cannot be seen by eye; 
however, since this is significantly smaller than the region below £ = 1 , it is not picked up when sampling linearly the 
parameter £. As we can see from the other two panels of figure [ 6 ] the viability priors do not have a strong impact on 
A and 77 in the region where the posterior of the model is peaked. 

From the top panels of figure [ 6 j we can notice that the different parameters specifying Horava gravity do not have 
sizable degeneracies between them over all the range explored (which spans several orders of magnitude). From the 
lower panel of the same figure we can notice that different data sets contribute differently to the parameter bounds. In 
particular we can see that PLC strongly constrains the A parameter, while preferring a bigger value for 77 and £. The 
addition of background probes pushes these two parameters closer to the ACDM limit of the theory. Noticeably the 
addition of CMB lensing strongly degrades the bounds on A while being consistent with PLC+BG for the other two 
parameters. This behavior is expected considering the known tension between the Planck 2013 and the Planck 2015 
data and the CMB lensing power spectrum as reconstructed from the CMB trispectrum [77) [92l [93] ■ These results are 
confirmed by the marginalized bounds on the H3 parameters reported in table [l] In particular the 99.7% C.L. bounds 
on the £ and 77 parameters weakly depend on the data set used and in particular for the 77 parameter no 99.7% C.L. 
bounds are found. As discussed in section 3.1 the “bare” cosmological constant A has been replaced in this analysis 


by Ll° DE , and the latter has been included in table [l] The most relevant result that can be drawn from table [l] is that 
low-energy Horava modifications of gravity are severely constrained, with the corresponding parameters bounded to be 
orders of magnitude below unity. In particular, we find that cosmological data have a strong constraining power on A. 
Here we summarize the bounds, at 99.7% C.L, that we get from the combination of all data sets: 


£ -1 = -o-oi±o:o2, 

l°gio(A — 1) < -4.31, 

f l° DE = 0.69 ± 0.02 (49) 


Let us notice that the viability priors give also an upper bound on £, i.e. £ < 1. All the marginalized constraints 
on standard and derived cosmological parameters for the H3 case are shown in appendix |B| 

In table[l]we report also the marginal bounds on the PPN parameters {an, 0 : 2 } and G cosmo /G n — 1. From the first 
two, we can notice the extreme complementarity of cosmological and solar system experiments in constraining Horava 
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Figure 7: Results of our analysis of the H2 case. Panel (a): The marginalized likelihood of log 10 (A — 1); Panel (b): 
The marginalized joint likelihood of log 10 (A — 1) and log 10 ? 7 . The darker and lighter shades correspond respectively 
to the 68 % C.L. and the 95% C.L.. The theoretical viability condition is roughly (A — 1) > 0(2r]/9)] Panel (c): The 
marginalized likelihood of log 10 ? 7 . In all panels different colours correspond to different data set combinations as shown 
in legend. The dark grey shade corresponds to the marginalized region of parameter space excluded by viability priors. 
The arbitrary normalization of the posterior is chosen to better display the effect of the viability priors. 


gravity. The cosmological observations lead to an upper bound on aq that is aq < 0.220, while PPN bounds on this 
parameter are three orders of magnitude stronger. Cosmological observations are, however, sensitive to £ and weakly 
sensitive to 77 , while solar system probes constrain just a degenerate combination of the two. On the other hand the 
cosmological bounds on the parameter a .2 are worst than the PPN ones by several orders of magnitude. That is due to 
the fact that in the limit in which A is constrained to be smaller than the other parameters by cosmological data, the 
a 2 parameter goes to infinity as it is clear from its definition (|35|) . If we compare the cosmological constraint on A with 


the one that is derived from a 2 in eq. (36), we see that our bounds are compatible with the solar system constraints. 
Finally, we can see that the bound on G cosmo /Gjv — 1 < 0.028 (99.7% C.L., all combined) is improved by one order of 
magnitude with respect to previous results HU. 


4.3 H2 case: results 

The second case we consider, is a sub-case of the previous one obtained by restricting to the plane of the parameter 
space corresponding to 77 = 2£ — 2. For this choice, the solar system constraints are automatically evaded as shown by 


studies of the PPN limit of Horava gravity (see section 2.4). We refer to this as the H2 case. The free parameters of the 


theory are now { 77 , A} and, as discussed in the previous section, we sample the parameter space of the logarithm of p 
and A — 1 since both these quantities are constrained by the viability priors to be positive. Unlike for the H3 case, where 
the viability priors had a strong influence only on the parameter £, in the H2 case they have a strong influence on the 
marginalized posterior for both 77 and A, as one can see in figure [Tj In particular one can notice in panel (b) of figure [7] 
that the viable region is a triangle corresponding roughly to the condition (A — 1) > 0(2ij/9). This triangle shape of 
the marginalized joint posterior of the model parameters has a strong influence on the marginalized ID posterior of 
the two parameters. In particular, as we can see in panel (a) of figure [7], the low tail of A is cut by the viability priors 
and, in panel (c), the posterior of 77 becomes markedly non-gaussian. Apart from the degeneracy induced by this prior 
cut, no other degeneracy between the Horava gravity parameters is present. 

In panel (a) of figure [7j we can see that different data sets contribute differently to the cosmological bounds. In 
particular, the PLC, PLC+BG and PLC+WiggleZ data sets are pushing the posterior of A to smaller values. The 
constraints become slightly weaker when considering CMB lensing for the same reasons explained in the previous 
section and it dominates the bounds coming from the total data set combination. Because of the degeneracy induced 
by the viability priors, in the H2 case we are able to set bounds also on 77 . These are shown in table [2] In this table we 
can notice that the 99.7% C.L. confidence bounds on A are comparable to the previous ones and the bounds on 77 are 
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H2 case 

Parameters 

Prior 

PLC 

PLC+BG 

PLC±CMBL 

PLC±WiggleZ 

all combined 

logi 0 (A - 1) 
logio V 

^DE 

[-10,-2] 

[-10,-2] 

derived 

< -4.62 

< -4.68 

0-68±g;“ 

< -4.62 

< -4.58 

0.69 ±0.02 

< -4.38 

< -4.40 

0.69 ± 0.04 

< -4.59 

< -4.73 

0 66 +0 05 
u - uu -0.06 

< -4.39 

< -4.51 

0.69 ±0.02 

^cosmo/^iV 1 

derived 

< 3.6 x 10 _b 

< 3.6 x 10~ b 

< 6.3 x 10 _b 

< 3.9 x 10 _b 

< 6.1 x 10" 5 


Table 2: The 99.7% C.L. marginalized posterior bounds on the H2 case parameters and G cosmo /GN — 1. 


considerably stronger. Indeed the bounds for all combination of data sets for the H2 case are: 

log 10 (A-l) < -4.39, 
log-i 0 (>/) < -4.51, 

Q° DE = 0.69 ± 0.02. (50) 

As in the previous case &P DE and all the standard derived cosmological parameters are reported in appendix |b| 

Noticeably in the H2 case, the bounds on G cosmo /GN —1(< 6.1 x 10 —5 , 99.7% C.L., all combined) are more stringent 
than in the H3 case. In particular this bound is several orders of magnitude stronger than the BBN bound. 

The H2 case studied in this section has been already investigated in ref. m Indeed, our H2 case is the khronometric 
model in ref. m when the LV in matter is switched off and a precise mapping between the different notations adopted 
is worked out. In this respect, a comparison between our results and the one in ref. m is not straightforward. We 
find an overall good agreement in the cosmological predictions by comparing the output of the code used by ]15i to 
EFTCAMB. The stronger bound on G cosmo /Gjv — 1 is then due to the different choices of cosmological parameters 
that we sample and the different priors that we impose on them. 


5 Conclusion 


In this paper we have performed a thorough investigation of the cosmology of Hofava gravity, which is a Lorentz 
violating theory proposed as a candidate for quantum gravity muz!. The emergence of Lorentz violations at all 
scales, makes the theory power-counting renormalizable at very high energies (~ 10 19 GeV/c 2 ). However, since at low 
energies we do not experience LV effects, we expect large scales tests to place important constraints on the theory. In 
particular, we have analysed Hofava gravity at cosmological scales, to see whether there is any room for LV at these 
energies. As we will summarize in details in the following, we have found that cosmological data severely constrain 
Hofava gravity. 

We have performed our analysis within the EFT framework for dark energy and modified gravity HOGS], which 
we have reviewed in section |2.2| We have focused on the dynamics of the background and linear perturbations, and 
considered the power-counting renormalizable action for Hofava gravity which includes all the quadratic operators with 
up to sixth spatial derivatives m- For this action, we have worked out a complete mapping into the EFT language, 
in section 2.3 finding that its low-energy part is completely mapped in the most commonly used EFT action ([7|. 
While the high-energy part requires the inclusion of additional operators, that we have identified in appendix [A] After 
working out the full mapping, when exploring the cosmology and corresponding observational bounds, in sections [3] 
and [4] we have restricted to the low-energy part of the action which is sufficient for a first exploration of the large scale 
phenomenology of the theory. 

For our analysis we have made use of the powerful EFTCAMB/EFTCosmoMC package [ID [421 S3]- While this 
package was made publicly available by some of the authors at http: //wwwhome. lorentz. leidenuniv.nl/~hu/codes/, 
for our analysis we have used an updated version which fully implements Hofava gravity. By the latter we mean that, for 
the first time, we have implemented in EFTCAMB a full mapping case, solving the Hofava equations for the background, 
instead of using a designer approach to the expansion history. We have included a detailed discussion of this procedure 
in section 3.1 After solving the background and before proceeding to the evolution of the perturbations, EFTCAMB 


runs a check on the viability of the selected theory, enforcing some stability requirements such as the absence of ghosts 
and gradient instabilities. The latter conditions become viability priors when using EFTCosmoMC to constrain the 
parameters of the theory by means of cosmological data. Finally, we have proceeded to evolve linear scalar perturbations 
with the general EFT equations [H| , specializing their coefficients to the corresponding expressions in the Hofava case 
through the mapping worked out in section [2.3[ 

Let us notice that besides the physical stability requirements, Lorentz violations can also be constrained via BBN and 
solar system tests, as discussed at length in section [2~fj In our analysis we have not imposed the BBN experimental 
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bounds a priori, rather we have compared them to our finding for cosmological constraints. Nevertheless, we have 
investigated two cases: a first one, H3, where the low-energy Horava gravity parameters {A, £, 77 } were allowed to vary 
freely; the second case, H2, where we enforced a relationship between the parameters that allows the theory to evade 
PPN constraints, reducing the number of free parameters to two, {A, 77 }. 

In section [3j we have studied in details the cosmology of Horava gravity. At background level, we have found a 
constant rescaling of the Hubble rate which reflects in the behaviour of the density parameters. Indeed, as shown 
in figure [T the fractional matter density exceeds unity at all times and the fractional density of the effective dark 
energy (44) correspondingly becomes negative, so that the flatness condition is satisfied at all times. This behaviour 
signals that the modifications of gravity in this theory should be considered as a purely geometrical effect, rather than 
be interpreted in terms of a dark fluid. At the level of linear perturbations we have identified modifications in the 
ISW effect, the gravitational lensing, the rate of growth of structure and the B-modes spectrum, which translated 
into stringent bounds on the Horava parameters when we fit them to cosmological data. In order to facilitate a in 
depth visualization of characteristic features of Horava gravity at the level of perturbations, we have specialized to two 
choices of parameters, one for the H3 case and one for the H2 case, and have performed a thorough analysis of the 
dynamics of perturbations and the corresponding cosmological observables. In both cases, we have found a general 
enhancement of the growth of matter perturbations and the lensing potential. The first modifies the shape of the matter 
power spectrum, which we have found to be enhanced for 10 ~ 3 < k < 10 ” 1 h/Mpc; the second one impacts the CMB 
temperature power spectrum at high-f? and the CMB lensing power spectrum. In particular in the latter case the effect 
is of a general enhancement of power in the lensing potential auto spectrum. A modification of the lensing potential 
also alters the scalar perturbation induced component of the B-mode power spectrum, leading to an enhancement in 
that signal. On the other hand, we found an enhancement of the CMB temperature power spectrum at large angular 
scales and its cross-correlation with the lensing potential. This effect can be traced back to an enhancement of the 
ISW effect at late times. Finally, the tensor power spectrum is also modified as the speed of sound of the tensor modes 
depends on the parameters of the theory. In particular the tensor BB-power spectrum is shifted on the right with 
respect to the ACDM one for the H3 case as the tensor propagation is sub-luminal and on the left for H2 because the 
propagation is super-luminal. Noticeably, we have determined that in general the quasi-static approximation is not safe 
to describe the evolution of sub-horizon perturbations in Horava gravity if we want an accuracy better than 30%. Let 
us stress that while the direction and entity of these modifications is dependent on the specific choice of parameters, 
we have found a general enhancement of the growth rate, lensing potential, and an enhancement of the ISW effect for 
several choices of parameters that we have sampled in the region allowed by the viability priors. 

In section [4j we have moved on to perform a global fit of the two cases of low-energy Horava gravity, H3 and H2, to 
progressive combinations of cosmological data sets: the CMB temperature-temperature and lensing power spectra by 
Planck 2013, WMAP low-1? polarization spectra, the WiggleZ galaxy power spectrum, the local Hubble measurements 
and Supernovae data from SNLS, SDSS and HST and the BAO measurements from BOSS, SDSS and 6 dFGS. 

In the case of H3, we have set upper bounds on A and a lower bound on £, while for (which through eq. (45) 
replaces A) we found a mean value and errors that are close to the ones of the ACDM model. Specifically we obtained 
£ — 1 = —O.Ollg pi,, log 10 (A — 1) < —4.31 and f l° DE = 0.69 ± 0.02 at 99.7% C.L. for the combination of all the data sets 
considered. For all the other data set combinations see table [l] As a general result we have found that the values of 
the Horava gravity parameters are constrained to be close to their values in the GR limit. For both cases we reported 
the constraints on the standard cosmological parameters in tables [3][4] Moreover, for the H3 case we get an improved 
bound on G cosmo /Gjv — 1 < 0.028 ( 99.7% C.L., for the combination of all data sets) which outruns the BBN one. On 
the other hand PPN experiments are three orders of magnitude stronger in constraining the on parameter that we find 
to be a\ < 0.220 ( 99.7% C.L., for the combination of all data sets), while our cosmological bound on A is compatible 
with the one derived from solar system tests. 

For the H2 case, we were able to set upper bounds on log 10 (A — 1) < —4.39 and log 10 ( 77 ) < —4.51 and constraints 
on GP de = 0.69 ± 0.02 at 99.7% C.L. with all data sets. Noticeably for this model we get a quite stringent bound on 
Gcosmo/G n — 1 < 6.1 x ICC 5 at 99.7% C.L. by combining all the considered data sets. 

The full mapping of the low-energy limit of Horava gravity has been publicly released as part of an update of 
EFTCAMB/EFTCosmoMC. As part of future work, it would be certainly of interest to explore the phenomenology as¬ 
sociated to the high-energy part of the Horava gravity action to see whether additional operators can affect significantly 
linear perturbations. Future analysis could also include the study of LV in the dark matter sector. 
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A The L 4 and Lq Lagrangians 

The L 4 and Lq Lagrangians contain, respectively, all the operators up to fourth and sixth order spatial derivatives, 
which are compatible with the symmetry of Horava gravity and guarantee its power-counting renormalizability. The 
number of those operators is very large, but given that we are interested in an effective field theory description of linear 
scalar perturbations, only the ones which are second order in perturbations have to be considered in the action. The 
latter have been identified in ref. and they are given by suitably contracting the three-dimensional Ricci tensor and 
scalar, the acceleration a,;, and their spatial derivatives. In detail, the HE part of action ([ 3 ]) can be written as 
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(51) 


where V 2 = V'V, and V 4 = V*ViV J Vj, and the coefficients gi are running coupling constants of suitable dimensions. 
The first and second lines contain respectively the contributions from L4 and Lq. 

We expand now the above action up to second order in perturbations by considering that on a flat FRLW background 
the components of TZ and 7 Zij identically vanish. Then we map it into the language at the basis of the EFT formalism 
discussed in section H by using the relation |l9| ) for a, : . With these recipes, it is straightforward to show that the 
operators in action (511 up to second order can be written as 
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We notice immediately that the EFT action |7| is incomplete if one wants to treat the full version of Horava gravity 
(which includes the operators in action ©), then we need to add to it all the operators in eqs. ([52]). This suggests to 
extend the EFT action discussed in section [272] to a more general one by adding the following part: 


5 


eft, 2 = J d A xV^ 9 [Mr)(Sn) 2 + X 2 (T)Sn. l] 6n^+ X 3 (T)6TZ\/ 2 {a 2 5g 00 ) + X 4 (T)d l (a 2 g 00 )V 2 d\a 2 g 00 ) (53) 

+ A 5 (t)6KV 2 6K + X e(T)\/ t 6n jk \/ i SU : > k + Xri^d^Sg^^id^^dg 00 )) + X 8 {T)\7 2 6K\7 2 {a 2 5g 00 )] . 


In the Horava gravity case the EFT functions A,;’s reduce to the constant coefficients in eqs. (52), e.g. Ai = 3 iTOq/( 2 £— 3 ). 

Notice that the first two operators in the action (53) have already been considered in ref. )33j , while the remaining 
operators have been considered in refs. [48] EO]. However, in these latter works an explicit EFT action (in the form of 
the action (531) has not been emphasized as well as an explicit mapping between these operators and a specific theory 
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has not been worked out. In this respect our finding corresponds to new results. Finally, let us mention that although 
we wrote the operators in eqs. (52) in terms of 3D quantities, following the 3+1 formalism employed in Horava gravity, 
one can always express them by means of 4D quantities by using the Gauss-Codazzi relation ;62l. 

It would be of interest to implement the contributions of these new operators in the equations for the perturbations 
evolved in EFTCAMB, in order to investigate their cosmological effects. We expect that their contribution becomes 
more important as the cosmological scale becomes smaller. This is part of ongoing work [1)4] . 


B Cosmological Parameters 

In this appendix we report the 99.7% C.L. constraints on the derived cosmological parameters: H 2 the present 

day density parameter of baryons; h 2 the present day value of the cold dark matter density parameter; 1000 mc 
which measures the sound horizon at last scattering; r is the Thomson scattering optical depth due to reionization; 
n s the power law index of the scalar spectrum; ln( 10 lo A s ) the log power of the primordial curvature perturbation, 
H 0 (km/s/Mpc) the present time value of the Hubble rate, and the dark matter density parameter today. We 
found that the constraints on these parameters for H2 are the same as in ACDM as reported in the following table. 
The reason for this is that, in H2 case, the Horava gravity parameters are constrained to be very close to the GR limit 
so that the cosmological parameters fall back to their ACDM values. 



Bounds on cosmological parameters for ACDM and H2 cases 

Parameters 

PLC 

PLC+BG 

PLC+CMBL 

PLC+WiggleZ 

all combined 

n° b h 2 

0.02202 ± 0.008 

0.02216 + 0.0007 

0.02210 + 0.0008 

0.02182 ± 0.0008 

0 0221 9+ u - uuuy 

u.uzzio_o 0007 

n° c h 2 

0.120 + 0.008 

0.118 + 0.004 

0 11Q+ 0 - 007 

U.ii^_o.oo6 

0.123 + 0.008 

0 11Q+ 0 - 004 
U.liy_ 0> 003 

1000mc 

1.041 + 0.002 

1.041 + 0.002 

1.041 + 0.002 

1.041 + 0.002 

1.041 ± 0.002 

T 

0 f)89 +0 ' 047 
u.uo^_ 0 036 

0 OQ2 +0,040 
U.U^z_ 0 034 

0-089i°;^ 

0 OS^ 0 ' 037 
U.UoO_ 0 033 

0 QQO +0 037 
u.uyu_ 0 032 

n s 

0 Q^Q - ^ 0 - 021 

u.yo j _ 0 020 

0.963 + 0.015 

0-96118:2!? 

0.953 + 0.021 

0.963 + 0.015 

ln(10 lo A s ) 

q noo+0.085 
O.Uoo_o 069 

8 OQO+ 0 - 078 
O.U^U_ 0 069 

8 088+ 0 - 073 
O.UoO_ 0 064 

q no7+0.072 

O.UO / —0 063 

3-o87i8;8J! 

Ho 

67.2 ± 3 3 l 

68.0l?;8 

67.71!;! 

65.8!^ 

67.7lfj 

oo 

^ L m 

0 81 f)+ 0 - 052 

U.Olu_o 044 

0 SO^ 0,022 

U.OUO_o Q24 

0 80Q +0 ' 042 

U.OUy_ 0 Q39 

O qq7+0-051 

U.OO / _o Q46 

0 80Q +0 ' 022 
u.ouy_ 0 021 


Table 3: Mean values and 99.7% C.L. bounds on several cosmological parameters. Notice that these bounds do not 
change for the ACDM and H2 cases. 



Bounds on cosmological parameters for H3 case 

Parameters 

PLC 

PLC+BG 

PLC+CMBL 

PLC+WiggleZ 

all combined 

n° b h 2 

0 021 +°- UUi 

u.uz ±_ 0 ooi 

n n9i8+ u -° uu9 

U.Uzlo_ 0 ooiO 

O 091 Q+9-9010 

U.UzIo_o 0009 

0.021 + 0.001 

n 09i q+9.uuu8 

U.Uzio_ 0 q 009 

n° c h 2 

0 12+ 001 
U, 1 Z - 0.01 

0 11Q+ 0 - 005 
U.±±9_ 0 Q04 

0.118 + 0.007 

O 1 oq+0.008 

u.izo_ 0 Q10 

0 11 Q+°- 004 

1000 mc 

1.04318:881 

1 Q4^+ 0 - 005 
i ' u ^°_0.003 

1 Q4Q+0.006 

l-U z ±9_o.o04 

1-04318:881 

1-04418:881 

T 

o.o8l8;8I 

0.08 + 0.04 

0.08 + 0.04 

0.08±8;8| 

o-o8i8;8| 

n s 

0.96 + 0.03 

0.96 + 0.02 

0.96 + 0.02 

0.96 + 0.03 

0.96 + 0.02 

ln( 10 lo A s ) 

o QQ+0.09 

o.uy_ 0 .08 

O qq+ 0.09 
°- uy -0.08 

3.08l°;° o 8 7 

3.10 + 0.08 

q rj 7 + 0.08 
° >u ' -0.06 

Ho 

67 + 5 

67 + 2 

68 I 3 

65 I 4 

67 I 2 

n° 

m 

0 Q1+0.07 
u - ol -0.05 

0.31 + 0.02 

o.3oir 0 i 

0.33 + 0.06 

0.30 + 0.02 


Table 4: Mean values and 99.7% C.L. bounds on several cosmological parameters in H3 case. 
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